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PACS 73 . 23 . -b - Electronic transport in mesoscopic systems 

PACS 73 . 20 . Fz - Weak or Anderson localization 

PACS 73 . 40 . -c - Electronic transport in interface structures 

Abstract. - We have studied numerically the mesoscopic fluctuations of the conductance of 
a graphene strip (width W larger than length L), in an ensemble of samples with different re- 
alizations of the random electrostatic potential landscape. For strong disorder (potential fluc- 
tuations comparable to the hopping energy), the variance of the conductance approximates 
the value predicted by the Altshuler-Lee-Stone theory of universal conductance fluctuations, 
VarGucF = 0.12 (W/L)(2e 2 /h) 2 . For weaker disorder the variance is greatly enhanced if the 
potential is smooth on the scale of the atomic separation. There is no enhancement if the poten- 
tial varies on the atomic scale, indicating that the absence of backscattering on the honeycomb 
lattice is at the origin of the anomalously large fluctuations. 



Phase coherent diffusion in metals is accompanied by 
sample-to-sample fluctuations in the conductance of the 
order of the conductance quantum e /h, dependent on the 
shape of the conductor but independent of its size or of 
the disorder strength. This is the phenomenon of univer- 
sal conductance fluctuations (UCF) [1,2]. The universality 
does not extend to different transport regimes, in partic- 
ular the fluctuations become much smaller than the UCF 
value both in the ballistic regime of weak disorder and in 
\ the localized regime of strong disorder [3,4]. 

In this paper we report on the observation in a computer 
simulation of a transport regime with conductance fluctu- 
ations that are much larger than the UCF value. The 
anomalously large fluctuations appear in a tight-binding 
' model of a carbon monolayer, for a disorder potential that 
is smooth on the scale of the atomic separation and weak 
on the scale of the hopping energy. It is known that such 
a potential in a honeycomb lattice can deflect the elec- 
trons but cannot scatter them backwards [5,6]. The con- 
sequences for weak localization of the absence of backscat- 
tering have been studied theoretically [7-10] and exper- 
imentally [11, 12]. While conductance fluctuations as a 
function of magnetic field in a given sample have been ob- 
served experimentally [11,13,14], and analyzed in terms 
of the UCF theory, the anomaly found here in the sample- 



to-sample fluctuations has not been reported previously. 

We consider a disordered strip of graphene in the x — y 
plane, connected to ballistic leads at x — and x = L 
(see Fig. 1). The orientation of the honeycomb lattice is 
such that the edges at y = and y = W are in the zigzag 
configuration. We vary L and W at fixed aspect ratio 
(mostly taking a rather large ratio W/L — 3 to minimize 
the effects of edge states). The lattice Hamiltonian is 

h = 5>tf|i)oi +E[ c w r *) + i*x»i- (!) 

i,j i 

The hopping matrix element = — t if the orbitals \i) 
and \j) are nearest neigbors (with r s» 3eV), otherwise 
Tij = 0. The velocity v near the Dirac point equals 
v = ^\Z3ra/h s» 10 6 m/s, with a = 0.246 nm the lattice 
constant. 

The electrostatic potential contains a contribution C/g a tc 
from gate electrodes and a random contribution U lulv from 
impurities. The potential t/ ga t c vanishes in the leads x < 
and x > L and equals Uo in the strip < x < L. By 
varying Uq at fixed Fermi energy /too in the leads, we can 
vary the Fermi energy /to — n^, — Uq in the strip. We 
take /ioo = t/2 and compare the two cases Uq = Hoa => 
jLto = and Uq = => /j,q = /ioo. The first case is an 
undoped graphene strip, the second case is a heavily doped 
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Fig. 1: Top panel: Top view of the honeycomb lattice in a 
graphene strip, connecting two electrical contacts at a voltage 
difference (gray rectangles). The samples used in the simula- 
tion are much larger than the one shown here. Bottom panel: 
Potential profile along the strip, showing the fluctuations from 
the disorder. 



strip (but still at sufficiently small Fermi energy that the 
linearity of the dispersion relation holds reasonably well) . 

We generate a realization of the disorder potential by 
randomly choosing Ai mp lattice sites Ri, R2, . . ■ i?Ar imp out 
of the total number N to t — f V3 LW/a 2 of sites in the 
disordered strip, and by randomly choosing the potential 
amplitude U n at the n-th site in the interval (—5,6). We 
then smooth the potential over a range £ by convolution 
with a Gaussian, 



EAmp(r) 



AW 

E 



U n exp 



Rn\ 



2e 



(2) 



In the special case ( ««, A^ mp = N to t each of the lat- 
tice sites in the strip has a randomly fluctuating poten- 
tial. This is the Anderson model on a honeycomb lattice 
studied in Ref. [15]. We contrast this model of atomic- 
scale defects with the case £ = aV3 of a potential which 
is still short-ranged on the scale of the system size but 
which varies smoothly on the atomic scale. Such a poten- 
tial could be realized by screened charges in the substrate. 
(The Gaussian smoothing is chosen for computational con- 



venience, and we have checked that the results are not 
sensitive to the type of smoothing considered.) 

Wc quantify the disorder strength by the dimcnsionless 
correlator 



K = 



LW 1 



Ntot Ntot 



(hv 



)2 A/2 
1 Jv tot 



^^Wfill/kptrj)) (3) 



of the random impurity potential (with vanishing average, 
(fimp) = 0). This single number K is representative on 
length scales large compared to the correlation length £. 
For the model potential (2) we find (for W) 



lV3(8/T) 2 (N imp /N tot )K 2 , 

1 N imp jv tot 

-EE cx ? 



""P n=l i=l 
1 



Rn I 



2£ 2 



if £ < a, 



^7rV3(C/a) 2 if £> a. 



(4) 



(5) 



For large fio the correlator Kq determines the transport 
mean free path l iT in Born approximation [6,7], 



k F K 



2 if ^ > a, 
1 if £ < a, 



(6) 



where kp = \^\/hv is the Fermi wave vector in the strip 
(which should be » 1/Z tr for the Born approximation 
to hold). The factor-of-two increase in ? tr for smooth 
disorder is due to the absence of backscattering in the 
honeycomb lattice [5,6]. The corresponding "classical" 
conductivity (without quantum corrections) is given by 

Cclass = (2e 2 /h)k F hr- 

Wc calculate the transmission matrix t numerically by 
means of a recursive Green function algorithm. The 
conductance G then follows from the Landauer formula 
G = (2e 2 //i)Trii t . (The factor of two accounts for the 
spin degeneracy.) By repeating the calculation for some 
300-3000 realizations of the disorder potential, we obtain 
the average conductance (G) and the variance VarG = 
(G 2 ) — (G) 2 . Results are shown in Fig. 2 at the Dirac 
point (/io = 0) and in Fig. 3 at Hq = r/2. 

The Altshuler-Lee-Stonc theory of universal conduc- 
tance fluctuations (UCF) gives a variance [1,2,4] 



Var Gucf = C 



se 



W 
L ' 



if W^>L, 



(7) 



with G = (3/tt 3 )C(3) = 0.116 and ((x) the Ricmann zeta 
function. For atomic-scale disorder, the symmetry index 
(3 — 1 (orthogonal symmetry) and the degeneracy factor 
s = 2 (only spin degeneracy). For smooth disorder, one 
has (3 = 4 (symplectic symmetry) and s = 4 (both spin 
and valley degeneracy). In each case, the variance thus 
has the same value VarG UCF = C {W/L)(2e 2 /h) 2 . 

In Figs. 2b, 3b we see that the conductance fluctuations 
approach the UCF value for sufficiently strong disorder. 
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Fig. 2: Average and variance of the conductance as a func- 
tion of the strength of the disorder potential, quantified by 
the correlator (3). These plots are for the case that the dis- 
ordered strip is at the Dirac point (/xo = 0). The data points 
are for a smooth, short-range impurity potential (correlation 
length £ = av3), with different values of the impurity den- 
sity iVjmp/iVtot- Open symbols are for L — 40 a, filled symbols 
for L = 72 a (at fixed aspect ratio W/ L — 3) . The solid and 
dashed lines are the Anderson model of atomic scale disorder 
[£ = 0, Mmp = iVtot, L = 40 a (solid) and L = 72 a (dashed)]. 

This is by itself remarkable, since the Altshuler-Lee-Stone 
theory requires weak disorder, such that the conductiv- 
ity a = (G}L/W >• e 2 /h. Our numerical data for VarG 
only approaches Var Gucf when the disorder is so strong 
that a ~ e 2 /h. For weaker disorder, the conductance fluc- 
tuations first rise to a peak value VarGp Ca k well above 
Var Gucf, and then drop to zero upon entering the bal- 
listic regime. 

The increase of the conductance fluctuations above 
the UCF value does not happen for the Ander- 
son model of atomic-scale disorder (solid and dashed 
curves). For smooth disorder the enhancement factor 
Var Gpeak/Var Gucf increases with increasing Fermi en- 
ergy no — it is therefore not restricted to the vicinity of 
the Dirac point. The enhancement factor also increases 
with increasing £ (not shown), but at fixed £ it is insensi- 
tive to the system size (compare open and filled symbols 



Fig. 3: Same as Fig. 2, but now for the case that the disordered 
strip is away from the Dirac point (/io = r/2 = fioo)- 

in Figs. 2b, 3b). The anomalous enhancement does not, 
therefore, appear to be a finite-size effect. 

The transport mean free path (6) at /io = t/2 is hr = 
4\/3a/Ko (for smooth disorder), so ln/L w 0.1 at the peak 
of maximal conductance fluctuations in the largest system 
considered. We are therefore well outside of the ballistic 
regime, but the UCF value characteristic of diffusion is 
not reached until the mean free path has been reduced 
by another factor of ten. By comparing the data in Fig. 
3a for L/a = 40 and L/a — 72, we can conclude that 
the diffusive regime (with a scale invariant conductivity) 
is not reached until Kq > 10 for smooth disorder, while 
the diffusive regime is reached already for Kq > 0.1 for 
atomic-scale disorder. 1 

While the disappearance of the anomaly for atomic- 
scale disorder unambiguously indicates that the symplec- 
tic symmetry of the Dirac Hamiltonian is responsible for 
it, we have not been able to explain our simulations con- 
sistently in terms of existing transport theories for Dirac 
fermions [18-26]. Certain features of the data suggest a 

1 ~We have calculated the third and fourth cumulants, to search 
for deviations from a Gaussian conductance distribution. In the 
parameter range of Fig. 3 no significant deviations are obtained in 
the case of smooth disorder. We need atomic-scale disorder to obtain 
significantly non-Gaussian distributions at large disorder strengths. 
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partial explanation. 

First of all, at the Dirac point (/io = 0), the enhance- 
ment of the conductance fluctuations happens in the same 
range of disorder strengths as the enhancement of the con- 
ductivity above the ballistic value [16, 17] 

4 e 2 

^ballistic — 7~ ■ (8) 

7r n 

This increase of a was predicted by Titov [26] as a man- 
ifestation of resonant transmission of evanescent modes. 
We would expect such transmission resonances to enhance 
the mesoscopic fluctuations, but we would also expect the 
effect to diminish as the evanescent modes become prop- 
agating away from the Dirac point. Instead, the peak in 
VarG becomes larger with increasing fi , while the peak 
in a disappears. 

A second striking feature of the numerical data is that 
an increase of the impurity density Ni mp /N to t and a de- 
crease of the impurity potential S at fixed K has no sig- 
nificant effect on the conductance (compare the different 
open symbols in Figs. 2,3, which all lie approximately on 
a single curve). This signifies that the transition from 
the anomalously large fluctuations at weak disorder to the 
UCF value at stronger disorder is not related to the Born- 
Unitarity transition of Ref. [23] (which should appear at 
smaller Kq for smaller -/V; mp /./Vtot)- 

The percolation transition of Ref. [24] is more likely to 
be at the origin of the strong increase of the conductance 
fluctuations away from the Dirac point (where /cf£ ^ 1): 
in the regime 0.1 < Kq < 10 in between the ballistic 
and diffusive transport regimes. One would expect the 
presence or absence of a percolating trajectory to produce 
large sample-to-sample fluctuations in the conductance, 
which would increase both with increasing kp and with 
increasing £ — as observed in our simulations. This inter- 
pretation would imply that the conductance fluctuations 
result from variations in trajectories rather than fluctua- 
tions in phase shifts. 

To support this interpretation we compare in Fig. 4 
the variance VarG of the sample-to-sample fluctuations 
with the variance Var M G of the fluctuations obtained in 
a given sample upon varying the Fermi energy jiQ over a 
narrow interval. The former quantity contains contribu- 
tions both from variations in trajectories and variations 
in phase shifts, while in the latter quantity variations in 
phase shifts give the dominant contribution. To improve 
the numerical efficiency we took W/L = 1 for this com- 
parison. The results for Var G (filled symbols in Fig. 4b) 
are similar to those plotted in Fig. 3b for W/L — 3: A 
large enhancement appears of the sample-to-sample fluc- 
tuations above the UCF value. In contrast, the variance 
Var M G of the energy-dependent fluctuations (open sym- 
bols) does not show this enhancement, instead agreeing 
well with the UCF prediction [which for W/L = 1 equals 
VarGucF = 0.186 x (2e 2 /h) 2 ]. 

In the Altshulcr-Lcc-Stone theory of UCF one has 
Var G = Var^ G: Sample-to-sample fluctuations and fluc- 
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Fig. 4: Average and variance of the conductance away from the 
Dirac point (fio — r/2) as a function of disorder strength, for 
a square sample (L = W = 121 a). The filled symbols in (b) 
give the variance Var G of the sample-to-sample fluctuations, 
while the open symbols give the variance Var M G of the energy- 
dependent fluctuations. (The latter quantity was calculated 
for a given sample by varying /io £ (0.44 t, 0.5 r) and then 
averaging the resulting variance over 400 samples.) The solid 
curves in (a) and (b) represent, respectively, (G) and VarG in 
the Anderson model of atomic disorder. (In that model there 
is no significant difference between VarG and Var M G.) 

tuations as a function of energy or magnetic field give 
the same variance. Our computer simulations imply that, 
remarkably enough, this socalled ergodicity of the meso- 
scopic fluctuations does not hold in graphene. An analyt- 
ical theory to explain this unexpected numerical result is 
still lacking. 
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